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Abstract 

We study two-colour QCD with two flavours of Wilson fermion at non-zero 



in 
o 

^ chemical potential fj, and zero temperature. We find evidence of two separate 

CN ■ transitions: an onset transition at \x k, m 7T /2 where quark number and energy 

!> . densities increase from zero; and a deconfmement transition at higher fi. The 

^-dependence of the number and energy densities and the diquark condensate 
indicate that a Fermi surface is formed and that BCS rather than Bose-Einstein 
condensation dominates at the quark mass considered here, especially beyond 
the deconfmement transition. 



1 Introduction 



QCD with gauge group SU(2) and non-zero chemical potential fi has a number of features that makes it 
an interesting object of study. Since the quarks and antiquarks live in equivalent representations of the 
colour group and can be related by an anti-unitary symmetry (the Pauli-Gursey symmetry), the fermion 
determinant remains real at non-zero \i. For an even number of flavours, it can also be shown to be 
positive, making standard Monte Carlo sampling possible. Thus two-colour QCD can be used as a 
laboratory for investigating gauge theories at high density and low temperature, a region of phase space 
Oh. which is still inaccessible to lattice simulations of real (three-colour) QCD. 

At n = 0, the Pauli-Gursey symmetry implies an exact symmetry between mesons and diquarks, 
which are the baryons of the theory. In particular, chiral multiplets will contain both mesons and baryons. 
For Nf = 2 for example, the pseudo-Goldstone multiplet consists of the pion isotriplet plus a scalar 
isoscalar diquark-antidiquark pair. The diquarks can be expected to condense when p > m n /2, forming 
a superfluid ground state. In this respect, the theory is radically different from real QCD, where no gauge 
invariant diquark condensate exist and the ground state at high density is superconducting. The nature of 
the superfluid ground state is however an interesting issue in its own right. 

When <C m p (p denoting any non-Goldstone boson) the /U-dependence of the theory can be 
analysed using chiral perturbation theory Q. In this case a transition to a Bose-Einstein condensate 
(BEC) of tightly bound diquarks occurs at \i = [i = m n /2, giving rise to a quark density n q and 
diquark condensate (qq), 




n q K /ttG" - Mo); (qq) oc \ — ) lim^(qq(fj,)) = const. (1) 

This has been confirmed in simulations with staggered fermions Cl|3]|4ll3. 

This picture is expected to break down when w m p , at which point the Goldstones are no 
longer distinguished hadrons, and the fermionic nature of the constituents comes into play. At this point 
one may expect BCS condensation of quarks near the Fermi surface to become the dominant mechanism. 
Assuming states within thickness A <C \i around the Fermi surface participate in the pairing, we obtain 
the following qualitative behaviour: 

n q oc fi 3 ; (qq) oc Ap 2 . (2) 
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The same rule of thumb predicts the quark energy density e q oc /A 

In the gluon sector, the differences between SU(2) and SU(3) are expected to be less important, 
and 2-colour QCD is a good setting for ab initio studies of gluodynamics in the presence of a background 
baryon density. Of particular interest is the issue of deconfinement at high density. Signals of deconfine- 
ment were observed in simulations with Wilson [6| and staggered fermions [7], where correlations were 
found between the Polyakov loop and chiral or baryonic observables. However, the phase structure has 
not been investigated in any further detail, and it remains unclear whether there is a confined phase with 
non-zero baryon number (as in QCD), or just a single phase transition. This will be investigated in the 
present paper. 

2 Lattice formulation 

The Nf = 2 fermion action is given by 

S = ^iM(m)^i + &M(/i)^ - J$x{CTs)v$t + J^iCl^u 0) 
where M(fj,) is the usual Wilson fermion matrix 

M xy (fi) = 6 xy - «£[(1 - lv )e^U v {x)8 y , x+i> + (1 + lu)e-^°Ui{y)5 y ^] . (4) 

The diquark source terms J, J serve a double purpose in lifting the low-lying eigenmodes in the super- 
fluid phase, thus making the simulation numerically tractable, and enabling us to study diquark conden- 
sation without any "partial quenching". The results will at the end be extrapolated to the physical limit 
J = J = 0. We will also introduce the rescaled source strength j = J/n. 

The Pauli-Giirsey symmetry is 

KM^K' 1 = M*(jm) with K = C^t 2 . (5) 

This symmetry implies that det M(ji) is real, but not necessarily positive. However, with the change of 
variables <j) = —^Cri, 4* = C~ l T2"4>2 we can rewrite the action as 

It is now straightforward to show that det M. is real and positive if J = J* lUE), and it is possible to 
simulate this action with a standard HMC algorithm ||9)- 

We have computed the quark number and energy densities, diquark condensate, gluon energy 
density and Polyakov loop on each trajectory. The quark number density is given by 

n i = Yl K (V^)(7o - l)e^U t {x)^i{x + t) +^(x)( 7o + ly-PUlix - t)ipi(x - t)) . (7) 

i 

The energy density is given by 

£<i = J2 K (^( a; )(7o - l)e»U t (x)^(x + t) - fc{x)(rf + l)e^uj(x - i)^{x - t)) - e° q , (8) 

i 

where the vacuum subtraction e° q is 

s° q = ^(NfN e - <V^)m=o) ■ (9) 

The diquark condensate is 

(qq) = ^ 1 K^ -^Kfa), (10) 

while the gluon energy density is given by e g = ^-(Tr(P t — P s )), where P t and P s are the timelike and 
spacelike plaquettes respectively. 
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3 Results 

We have simulated with (3 = 1.7, k = 0.178 on an 8 3 x 16 
lattice, corresponding to a = 0.22 fm, m^/nip = 0.92 [ 10 1 . We 
have used a diquark source aj = 0.04 with 0.3 < a/x < 1.0. 
The number of trajectories and their length are given in tabled 
We have also started simulations for selected values of /i with 
aj = 0.02, but those results will not be shown here. Configu- 
rations were saved every 4 trajectories of length 1 and every 8 
trajectories of length 0.5. 

In figure [2 we show results for the Polyakov loop, gluon 
and quark energy densities, quark number density and diquark 
condensate, as function of chemical potential. We see evidence 
of two distinct transitions: first, an onset transition around a/i = 
0.4 where the quark number and energy density increase from 
zero; and second, a deconfinement transition around a/i = 0.65 
where the Polyakov loop takes on a non-zero value. We also 
find that the chiral condensate starts to decrease from its vacuum 
value around afi = 0.4. 

appears strikingly different 
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The behaviour of n q and 



Table 1: Run parameters for the HMC sim- 
ulations. N is the number of trajectories, 
while £ is the average trajectory length, 
from that of (0, and is suggestive of the formation of a Fermi surface and a /i-dependence qualitatively 
similar to that of © [9]. On the coarse lattices used here, the continuum behaviour may however be 
obscured by lattice artefacts. To compensate for this, we divide our results with the expected behaviour 
of free, massless Wilson fermions on an Nf x N t lattice, 



where 



WfN, 
N*N t 



E 



i sin ko [£V cos ki — 4] 
[4-^cos^] 2 + E,sm 2 ^ 




S(«o + |) 



i/j,. 



0, 

1,2,3. 



(11) 



(12) 



and N s = 8, N t = 16 in our case. The result of this is shown in figure 12 The number density shows 
good BCS scaling already from a/i 0.5, while no such agreement is shown for the energy density until 
the deconfinement transition. The fact that n q exceeds the free field value by a factor of two is consistent 
with a degenerate system with kp > Ep, suggesting a non-zero binding energy. Strikingly, the gluon 
energy density exhibits very good scaling with /x 4 over the whole range of \x considered. 

Also in fig. |3 we show the diquark condensate divided by /j 2 , which should be constant in the 
case of BCS condensation. Such scaling is indeed observed in the deconfined phase, for a/i > 0.65. 
Below this point the behaviour appears to be in qualitative agreement with the x?T expression (Q, but 
this requires further study before any firm conclusion can be drawn. 

Finally, in figure |3] we show the static quark potential for selected values of \i. We see clear 
evidence of screening for a\x > 0.4. However, for a/i > 0.7 a new pattern appears to emerge, in that 
the short-distance potential is strongly modified (and suppressed) while the long-distance screening, if 
anything, is reversed. This is most clearly seen in the right-hand plot of fig. |3] where the deviation of 
the potential from that at // = j = is shown. The source of this behaviour is subject of ongoing 
investigation. 

4 Discussion 

We have studied the /i-dependence of local observables and the static quark potential for 2-colour QCD 
with a quark mass well beyond the expected range of validity for chiral perturbation theory. We find clear 
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Fig. 1: Results for bosonic (left) and fermionic (right) quantities, as a function of chemical potential /i, for aj — 0.04. The 
open triangles are preliminary results for the diquark condensate with aj — 0.02. 
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Fig. 2: Left: Quark number density n 9 and energy density e q , divided by their values for free massless fermions on an 8 3 x 16 
lattice. Right: The diquark condensate (qq) divided by the BCS formula 0, together with the average Polyakov loop (L). 
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Fig. 3: The static quark potential for different values of the chemical potential fi. The right-hand plot shows the potential 
normalised to its value at jit = j = 0. 
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evidence of two separate (onset and deconfining) transitions, and strong indications that a Fermi surface 
is formed, leading to the formation of a BCS condensate at high densities. In this sense, the system is 
more QCD-like than 2-colour QCD at lighter quark masses, where BEC dominates. 

The simulations have been carried out with a fixed, non-zero diquark source j. We are currently 
performing simulations with a second, smaller diquark source to enable an extrapolation to j = 0. This 
will lay the grounds for drawing more definite conclusions regarding the nature of the system in the 
different regions. 

We are also studying the electric and magnetic gluon propagator, which may provide useful input 
into the gap equation at high densities, and are planning to study the hadron spectrum, focusing in 
particular on the fate of the vector meson, which in a previous study [ 6 1 was found to become lighter at 
higher density. 
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